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1 .  INTERACTION 

The  interaction  of  shock  waves  with  boundary  layers 
has  long  interested  researchers  in  fluid  dynamics,  because  it 
is  a  phenomenon  of  considerable  practical  importance  in  such 
problems  as  flows  over  wings,  control  surfaces,  through  engine 
intakes  or  across  transonic  cascades.  However,  despite  the 
extensive  work  carried  out  in  the  field,  there  remain  many 
fundamental  unanswered  questions  for  the  three  dimensional  case. 

In  particular,  the  interaction  of  a  boundary  layer 
with  a  perpendicular  oblique  (skewed)  shock  wave  has  attracted 
much  interest  both  because  of  its  geometrical  simplicity  and 
of  its  practical  importance.  For  practical  applications,  the 
incoming  boundary  layer  is  turbulent  and  this  is  why  most  of 
the  research  has  been  done  on  this  type  of  interaction 
(Refs.  1-8). 

A  series  of  studies  was  conducted  at  Princeton  Univer¬ 
sity  on  this  topic  with  the  sponsorship  of  AFOSR  and  other 
agencies  (Army  Research  Office,  Naval  Air  System  Command,  NASA). 
They  led  to  a  better  understanding  of  the  parameters  of  impor¬ 
tance  which  determine  the  scale  and  characteristics  of  this 
interactive  flow  field  (Ref.  9).  The  detailed  flow  field 
measurements  performed  by  Oskam  (Ref.  6)  also  serve  now  as  test 
data  for  comparison  with  numerical  solutions  using  full  com¬ 
pressible  Navier  Stokes  equation  solvers  (Ref.  10). 

Nevertheless,  the  fundamental  mechanism  of  the  inter¬ 
action  is  still  unknown  and  it  appears  relevant  to  go  back  to 
the  study  of  the  more  simple  laminar  interaction  case,  as  was 
done  for  the  two  dimensional  interaction.  Moreover,  the  testing 
of  the  computer  codes  solving  the  Navier-Stokes  equations  can 
more  suitably  be  done  in  the  laminar  than  in  the  turbulenct  case 
since  there  is  no  uncertainty  brought  out  by  the  turbulence  level. 
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Therefore,  a  program  of  studies  was  established, 
aiming  at  determining  the  basic  characteristics  of  the  laminar 
boundary  layer/skewed  shock  wave  interaction  in  order  to  dis¬ 
cover  its  key  parameters.  Surface  flow  visualizations  and 
surface  pressure  measurements  were  made  over  a  wide  matrix  of 
flow  parameters  (Reynolds  number,  boundary  layer  thickness, 
shock  strength).  The  results  of  the  tests  are  described  in 
detail  in  section  3.  In  particular,  the  variation  of  upstream 
and  downstream  influences  as  a  function  of  the  flow  quantities 
was  examined  and  compared  with  their  two  dimensional  or  turbu¬ 
lent  counterparts. 

Simultaneously,  theoretical  investigations  were  under¬ 
taken.  Integral  methods  have  proved  valuable  for  the  prediction 
of  2D  laminar  boundary  layer/shock  wave  interactions  (Ref.  11). 
Some  attention  was  given  to  the  possibility  of  extending  such  a 
method  to  the  computation  of  a  quasi  three  dimensional  plane 
flow  (see  section  4).  Recently,  various  numerical  techniques 
were  proposed  to  solve  the  full  compressible  Navier  Stokes  equa¬ 
tions  (Ref.  12-14)  among  which  is  the  implicit  corrected  vis¬ 
cosity  technique  proposed  at  the  VKI.  A  preliminary  investiga¬ 
tion  of  this  method  was  performed  and  the  results  are  presented 
i  n  section  4 . 

The  possibility  of  performing  Laser  Doppler  Veloci- 
metry  measurements  is  treated  in  section  5,  section  6  presents 
the  conclusions  of  this  report. 
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2 .  FACILITY  AND  INSTRUMENTATION 

2 . 1  Wi nd  tunnel  facility 

The  experiments  were  carried  out  in  the  supersonic 
wind  tunnel  S-l  at  the  von  Karman  Institute.  This  is  a  con¬ 
tinuous  closed  circuit  facility  of  the  Ackeret  type  with  a 
40  cm  x  40  cm  test  section.  The  range  of  stagnation  pressure 
is  from  0.1  to  0.3  bar.  The  Mach  number  2.25  nozzle  was  selec¬ 
ted  to  perform  the  experiments.  This  gives  unit  Reynolds  num¬ 
bers  in  the  ranges  106  to  3106/m.  The  models  were  at  near 
adiabatic  wall  temperature  for  all  tests. 

2 . 2  Test  model  and  instrumentation 

2.2.1  Genera ^  de^cr ip t_i o_n 

The  model  consisted  of  a  sharp  leading  edge  flat 
plate,  together  with  a  series  of  wedges  with  incidences  a  = 
4,6,8,10  and  15°  with  respect  to  the  free  stream.  It  is  shown 
in  figure  1.  The  wedges  were  15  cm  high  in  order  to  avoid  any 
interference  from  edge  effects  and  18  cm  long.  They  could  be 
fixed  on  the  plate  at  three  different  longitudinal  stations, 
i.e.,  with  their  apex  at  Xp  =  6,9  or  12  cm  from  the  flat  plate 
leading  edge.  This  provided  boundary  layer  thicknesses  ranging 
from  1.1  to  2.2  mm  with  unit  Reynolds  numbers  of  1.2  and  2.4  106/m 
Therefore,  the  Reynolds  number  based  on  5  was  varied  between 
1800  and  3600.  The  flat  plate  spanned  the  tunnel  and  was  60  cm 
long. 


The  flat  plate  was  instrumented  with  four  rows  of 
pressure  taps  parallel  to  the  free  stream  direction  (i.e.,  wind 
tunnel  axis)  at  lateral  positions  of  y  =  1,5,9  and  13  cm  from 
the  wedge  apex  location.  These  positions  were  selected  after 
examining  surface  flow  visualization  patterns  (see  chapter  3). 
The  diameter  of  the  pressure  taps  was  1  mm.  The  inside  diameter 
of  the  tubing  was  0.75  mm. 
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The  pressure  tappings  were  studied  using  12  port 
model  N0601/1P-12T  scanivalves  equipped  with  Validyne  DP15 
transducers.  The  scanivalves  were  automatically  controlled  by 
a  mi croprocessor  PET  Commodore  CBM  3032.  The  output  of  the 
transducers  was  also  treated  by  this  microprocessor  system. 

With  this  on-line  system  shown  in  figure  2.  tabulated  pressure 
distributions  could  be  obtained  immediately  following  a  test. 

2.2.2  Surface  fl  ow_v  i_si^al_i^ati  on£ 

Surface  flow  visualizations  were  performs  . i ng  an 
oil-graphite  mixture.  This  was  adapted  from  a  kero  e-graphite 

technique  used  at  Princeton  (Ref.  15).  To  do  this,  snorter 

flat  plate  was  built  in  order  to  fit  into  the  S-l  t  >1  injec¬ 
tion  mechanism.  Its  fore  part  was  coated  with  the  .  _raphite 
mixture  prior  to  the  test  and  it  was  then  injected  into  the 
tunnel.  The  oil  flowed,  leaving  on  the  surface  graphite  par¬ 
ticles  which  produced  the  visualization.  The  pattern  was 
recorded  by  lifting  it  off  using  transparent  tape  and  then 
sticking  it  on  white  paper.  On  figure  3  is  shown  a  photograph 
of  the  model  after  a  test. 

2.2.3  Sh^ck,  wa v_e_posH2°n 

As  mentioned  by  Dolling  et  al .  (Ref.  9),  it  is  impor¬ 
tant  to  have  an  accurate  measure  of  the  shock  angle  8  with  res¬ 
pect  to  the  free  stream.  This  is  not  very  different  from  the 
theoretical  oblique  shock  wave  but  slight  errors  in  the  value 
of  6  may  lead  to  significant  inaccuracies  in  shock  wave  position 
at  large  lateral  distances.  The  value  of  S  is  approximately 
constant  spanwise  expect  near  the  wedge  apex,  which  ensures  that 
the  shock  wave  is  effectively  plane  spanwise. 

The  position  of  the  shock  wave  was  measured  using  a 
pitot  probe  (Ref.  16)  and  locating  the  pressure  rise  through  the 
shock  wave.  A  lateral  traverse  mechanism  was  mounted  on  a  sting 


3.  EXPERIMENTAL  PROGRAM 


3 . 1  Surface  flow  visualizations 

Surface  flow  visualizations  were  obtained  for  the 
following  experimental  conditions  :  a  =  4°,  6°  and  85; 

Xp  =  6,9  and  12  cm;  Re^  =  1.2  10^/m  and  2.4  10r'/m.  Typical 
examples  are  shown  in  figures  4  and  5.  Only  the  flow  upstream 
of  the  separation  line  was  visualized  as  it  proved  to  give  bad 
results  to  introduce  oil  near  the  fin  to  visualize  the  flow 
behind  the  separation  line.  Additionally,  visualizations  were 
performed  at  a  =  10°  and  15°  in  order  to  check  the  laminarity 
of  the  interaction  under  critical  conditions.  This  will  be 
treated  later  in  section  3.5. 

Several  observations  can  be  made  when  examining  these 
visualizations.  The  first  and  probably  most  striking  one  is 
that  a  separation  line  can  be  observed  even  for  the  weakest 
wedge  incidence  (4°),  i.e.,  for  an  inviscid  pressure  ratio  of 
only  1.27.  This,  according  to  the  terminology  of  Peake  et  al. 
(Ref.  17)  is  a  global  separation  line  since  its  emanates  from  a 
saddle  point  .  This  global  character  is  further  proved  by  the 
absence  of  oil  behind  the  separation  line.  The  extent  of  the 
separated  region  is  considerable.  When  this  is  compared  with 
the  two  dimensional  wedge  angle  needed  to  induce  incipient 
separation  as  computed  using  the  formula  by  Riethmuller  et  al. 
(see  Ref .  18 ) . 

3  7\/M 

QO 

a  .  =  — - 

1  Re  0.22 
XF 


*  There  must  exist  a  saddle  point  because  locally,  the  wedge 
being  blunt,  there  is  a  nodal  point  of  attachment  on  it  and 
because 
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one  realizes  that,  for  the  experimental  conditions  used,  this 
varies  between  3.5  to  4.7°  depending  on  Xp  and  Reu .  Therefore, 
for  a  two  dimensional  interaction  yielding  the  same  inviscid 
pressure  ratio,  the  flow  would  be  attached  (for  small  ReXp)  or 
would  exhibit  a  small  amount  of  separation  (for  large  ReXp) 
whereas  in  this  case,  an  extended  separation  is  observed  in  all 
cases . 


When  comparing  the  shape  of  the  separation  lines  at 
different  wedge  incidences,  one  observes  that  the  Reynolds 
number  effect,  if  any,  is  small.  The  effect  of  shock  strength 
on  the  other  hand,  is  measurable  but  most  of  it  is  due  to  the 
outward  displacement  of  the  shock  wave  when  the  wedge  angle  is 
increased.  This  is  in  contrast  with  two  dimensional  interactions 
when  the  separation  length  is  strorgly  depending  on  shock  strength 
(Refs.  19,20).  For  reasons  related  to  model  construction,  the 
wedges  were  asymmetrical,  i.e.,  the  incidence  and  the  angle  of 
the  other  surface  with  respect  to  the  free  stream  were  not  the 
same.  This  has  no  effect  on  the  inviscid  flow  field  since  the 
two  flows  are  independent  of  one  another.  On  the  contrary,  it 
proved  to  have  an  effect  on  the  viscous  flow.  Indeed,  when  one 
observes  the  flow  visualization  for  a  small  wedge  incidence, 
the  separation  line  is  attached  to  the  wedge  apex  (Fig.  4). 

For  larger  incidences,  the  separation  line  is  ahead  of  the  wedge 
and  its  most  forward  position  does  not  lie  on  a  longitudinal 
axis  passing  through  the  wedge  apex,  due  to  the  asymmetry  of  the 
wedge.  This  seems  similar  to  the  phenomenon  of  shock  detachment 
in  inviscid  supersonic  flows  and  would  be  worth  theoretical 
investigation.  Such  a  phenomenon  can  also  be  perceived  in 
figure  6  taken  from  Korkcgi  (Ref.  2). 

In  three  dimensional  flow  fields,  there  are  some 
particular  quasi  two  dimensional  ones  which  are  labelled  coni¬ 
cally  or  cy 1 i nd r i c a  1 1 y  symmetric.  For  example,  the  three  dimen¬ 
sional  hypersonic  corner  flow  (Ref.  21)  is  of  the  conical  type 
(Fig.  7).  In  the  present  case,  the  separation  line  is  linear 
near  the  wedge  apex  (except  when  detached  from  it),  then 
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becomes  curved.  The  question  can  be  posed  whether  the  inter¬ 
action  reaches  cylindrical  symmetry  far  from  the  wedge  apex, 
i.e.,  whether  all  variables  become  constant  along  lines  parallel 
to  the  shock  wave.  This  question  cannot  be  answered  by  the  flow 
visualization  since  one  cannot  observe  any  region  where  the 
separation  line  runs  parallel  to  the  inviscid  shock  wave  but 
this  could  be  due  to  an  insufficient  size  of  the  model.  A 
theoretical  discussion  of  the  issue  will  be  presented  in 
chapter  4. 


By  varying  the  wedge  height,  it  was  verified  that  for 
the  largest  height,  the  flow  could  be  supposed  of  the  semi¬ 
infinite  type  in  the  region  of  interest  for  the  measurements. 

3 . 2  Experimental  determination  of  the 
shock  wave  position 


The  first  fact  that  appeared  during  these  measurements 
was  the  importance  of  the  height  of  the  fin,  expecially  at  large 
lateral  distances.  The  existence  of  possible  edge  effects  can 
be  estimated  by  a  simple  theory  which  is  explained  in  the 
following.  Consider  the  experimental  configuration  shown  in 
figure  8.  The  Mach  cone  originating  from  the  tip  of  the  wedge 
(labelled  A  in  the  figure)  intersects  the  shock  wave  plane  along 
a  straight  line  AB.  Since  the  flow  is  essentially  inviscid  in 
that  part  of  the  flow  f’eld,  it  can  be  stated  that  the  shock 
wave  is  plane  below  line  AB.  Now,  the  condition  for  studying 
a  semi -i nf i ni te  interaction  is  that  the  shock  wave  is  effec¬ 
tively  plane  when  submitted  to  interaction  with  the  boundary 
layer.  Since  line  AB  is  going  down  as  distance  is  increased 
away  from  the  wedge  tip,  there  will  be  a  limiting  distance  at 
which  the  shock  wave  is  no  longer  plane.  From  reference  22  it 
is  known  that  the  height  of  the  interaction  region  is  approxi¬ 
mately  s,  thickness  of  the  incoming  boundary  layer.  Since  this 
is  only  a  few  millimeters,  it  will  be  neglected  with  respect  to 
the  wedge  height  and  the  limit  distance  will  be  calculated  as 
the  one  which  corresponds  to  the  intersection  of  line  AB  with 


the  flat  plate.  This  can  easily  be  calculated  using  the  geomet¬ 
rical  relations  layed  out  in  figure  8.  Once  the  angle  y  of 
inclination  of  line  AB  is  known  {see  Fig.  8),  the  position  of 
its  intersection  with  the  flat  plate  is  easily  computed,  in 
particular  its  lateral  distance  from  the  wedge  root  R.  For 
h  =  15  cm,  Mx  =  2.21,  and  a  =  4,6  and  8°,  the  values  of  this 
lateral  distance  are  31.8  cm,  27  cm  and  24  cm.  These  are  well 
beyond  the  largest  distance  considered  in  the  pressure  measure¬ 
ments  (y  =  13  cm).  Therefore,  the  interaction  may  be  considered 
of  the  semi - i nf i ni te  type  for  all  test  conditions. 

A  first  series  of  experiments  was  conducted  by  explor¬ 
ing  the  flow  with  a  pitot  probe  at  different  heights  over  the 
flat  plate,  with  a  wedge  5  cm  high.  The  location  of  the  sharp 
pressure  rise  associated  with  the  shock  wave  never  reached  an 
asymtotic  value  as  the  height  over  the  flat  plate  was  increased, 
proving  that  the  interaction  was  not  of  the  semi -i nf i ni te  type, 
as  could  be  predicted  from  the  previous  analysis.  Subsequently, 
the  experiment  was  repeated  with  a  wedge  15  cm  high.  Then  an 
asymptotic  value  was  noticed  for  all  wedge  incidences  and  all 
Reynolds  numbers.  In  order  to  verify  that  these  values  were 
indeed  associated  with  the  inviscid  shock  wave,  the  experiment 
was  redone,  in  which  just  the  wedge  was  placed  in  the  tunnel, 
without  the  flat  plate  thus  avoiding  any  interference  from  the 
shock/boundary  layer  interaction.  Values  for  the  shock  locations 
in  agreement  with  the  previous  "asymptotic"  values  were  recorded. 
Nevertheless,  some  problems  remained.  First,  the  shock  wave 
angle  as  determined  by  the  measured  position,  constant  in  the 
spanwise  direction  except  near  the  wedge  apex,  was  greater  than 
the  theoretical  value  by  as  much  as  2°  for  the  lowest  Reynolds 
number  whereas  the  difference  was  never  greater  than  0.30  in 
Dolling's  experiments  at  Mach  3  (Ref.  9),  conducted,  however, 
at  much  higher  Reynolds  number.  Second,  there  was  a  marked 
effect  of  unit  Reynolds  number,  the  measured  shock  wave  angle 
being  reduced  by  1°  when  changing  Rey  from  1.2  106/m  to 
2.1  10D/m.  Third,  the  pressure  distributions  recorded  with  the 
pitot  probe  moving  axially  exhibited  some  strange  behaviour. 
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Instead  of  measuring  a  steplike  distribution,  the  pressure 
decreased  immediately  after  the  shock  wave  when  moving  downstrea 
This  could  not  be  attributed  to  edge  effects.  This  sheds  doubts 
about  the  sense  of  the  measure  and  raises  the  issue  of  the 
existence  of  a  systematic  error  due  to  the  measurement  technique 

The  results  being  inconclusive,  they  will  not  be  used. 
A  comparison  with  optical  visualization  could  help  clarify 
this  issue. 

3 . 3  Surface  pressure  measurements 

Surface  pressure  measurements  were  obtained  for  the 
same  experimental  conditions,  i.e.,  a  =  4°  ,6°  and  83  ; 

Xp  =  6,9  and  12  cm;  Rey  =  1.2  10°  and  2.4  105/m.  Pressure  dis¬ 
tributions  were  recorded  at  lateral  distances  y  =  1,5,9  and 

13  cm.  Typical  examples  when  the  pressure  is  non  dimensiona- 

lized  by  stagnation  pressure  ahead  of  the  shock,  P0,  are  rep¬ 
resented  in  figures  9  to  11.  Each  figure  shows  the  pressure 
distribution  at  differens  y's  for  given  a ,  Xp  and  Reu  .  Pressure 
distributions  are  shown  as  a  function  of  the  distance  x  from  the 
flat  plate  leading  edge.  Also  indicated  on  the  figures  are  the 
locations  of  the  theoretical  shogk  wave  and  of  the  separation 
line  as  observed  from  the  visualizations  together  with  the 
inviscid  pressure  values  upstream  and  downstream  of  the  shock 
wave.  Clearly,  from  the  measured  distributions,  it  can  be 

stated  that  the  pressure  tappings  did  not  cover  the  entire 

interaction  region.  In  order  to  obtain  complete  information 
about  the  flow,  new  tappings  should  be  designed  and  the  tests 
should  be  redone.  This  is  the  object  of  the  research  proposal 
submitted  as  a  continuation  of  grant  AFOSR  82-0051. 

By  examining  the  pressure  distribution,  some  general 
observations  can  be  made.  For  all  shock  strengths  and  all 
lateral  distances,  the  pressure  distribution  can  be  divided  in 
two  parts.  Upstream  of  the  inviscid  shock  wave,  there  is  a 
pressure  rise  constant  in  magnitude  whenever  the  lateral 


distance  or  Reynolds  number  are,  and  varying  little  with  shock 
strength.  This  is  followed,  downstream  of  the  inviscid  shock 
by  a  rather  sharp  pressure  rise.  The  curve  overshoots  the 
inviscid  pressure  value  as  calculated  from  oblique  shock  theory, 
then  relaxes  to  the  inviscid  value.  The  pressure  rise  in  the 
second  part  of  the  distribution  depends  thus  strongly  on  the 
shock  strength  value  since  the  overall  pressure  rise  must  be 
equal  to  the  inviscid  pressure  rise.  The  overshoot  in  the 
pressure  distribution  is  probably  due  to  an  attachment  pheno¬ 
menon.  A  similar  division  of  the  pressure  distributions  in 
two  parts  was  also  noticed  in  the  turbulent  case  (Ref.  6) 
when  separation  is  present. 

For  the  skewed  shock  wave  turbulent  boundary  layer 
interaction,  it  is  shown  in  figure  12  (taken  from  Ref.  6)  that 
for i nci dences  greater  than  93,  the  pressure  at  the  position  of 
the  inviscid  shock  wave  is  little  affected  by  a  change  in  shock 
strength  whereas  the  downstream  pressure  rise  changes  signifi¬ 
cantly.  Also,  when  comparing  pressure  distributions  in  isobar 
form  (also  from  Ref.  6)  taken  with  two  different  boundary  layer 
thicknesses  (Figs.  13,14),  it  can  also  be  noticed  that  the  pres¬ 
sure  at  the  location  of  the  shock  wave  does  not  vary  with 
lateral  distance  and  is  approximately  the  same  in  both  cases. 

In  pressure  distributions  encountered  in  the  blunt 
fin  induced  shock  wave  turbulent  boundary  layer  interaction. 
Dolling  et  al.  (Ref.  23)  also  noticed  the  division  of  the  pres¬ 
sure  distribution  in  two  parts,  the  first  of  which  consisted  in 
a  pressure  rise  whose  magnitude  varied  little  as  lateral  distance 
from  the  centerline  was  increased  or  as  the  fin  leading  edge 
geometry  was  altered.  Such  a  behaviour  thus  seems  to  be  a 
characteristic  feature  of  three  dimensional  interactions. 


In  the  present  work,  the  pressure  level  at  the  location 
of  the  shock  was  approximately  1.12  p2 •  This  is  of  the  same 
order  of  magnitude  as  the  pressure  at  separation  for  a  two  dimen¬ 
sional  interaction.  For  M  =  2  and  the  various  Reynolds  numbers 
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used,  the  latter  varies  between  1.12  and  1.17  pj  .  It  is  impor¬ 
tant  to  note,  however,  that  in  the  two  dimensional  case,  the 
pressure  level  at  separation  varies  strongly  with  Reynolds 
number  whereas  the  pressure  rise  in  the  three  dimensional 
separation  line  is  much  lower  than  the  latter  value  (see 
Figs.  9-11). 

Finally,  all  pressure  distributions  exhibit  a  kink 
near  the  inviscid  shock  wave  location.  This  kink  tends  to  be¬ 
come  a  trough  as  lateral  distance  Y  is  increased.  A  similar 
trend  in  also  observed  when  Reynolds  number  increases.  Such  a 
kink  is  generally  associated  with  the  existence  of  separation 
in  two  dimensional  interactions  (Ref.  20). 

A  useful  piece  of  information  which  can  be  syste¬ 
matically  extracted  from  the  pressure  distribution  is  the  extent 
of  upstream  influence.  In  order  to  render  more  precise  the 
determi nation  of  upstream  influence,  it  has  been  defined,  as 
indicated  in  figure  15,  the  straight  line  appearing  in  the 
drawing  being  the  tangent  to  the  pressure  distribution  at  the 
inflexion  point.  The  values  of  upstream  influence  in  functions 
of  Re,  a,  Xp  are  presented  in  Table  I.  From  the  examination  of 
the  dimensional  quantities,  several  observations  can  be  made. 

First,  there  is  hardly  any  effect  of  shock  strength 
on  the  amount  of  upstream  influence.  This  is  in  sharp  contrast 
with  two  dimensional  interactions  where  the  shock  strength  is 
the  leading  parameter  in  the  determination  of  upstream  influence. 
A  discussion  of  this  phenomenon  will  be  presented  in  chapter  3.6. 
On  the  contrary,  the  downstream  influence  depends  on  shock 
strength  :  it  decreases  as  shock  strength  increases.  This  is 
also  in  sharp  contrast  with  two  dimensional  interactions,  where 
downstream  influence  varies  in  the  same  sense  as  upstream  influ¬ 
ence  when  shock  strength  is  increased,  i.e.,  they  both  increase. 

A  possible  explanation  is  offered  for  the  decrease  in  downstream 
influence  as  shock  strength  is  increased  in  the  three  dimensional 
case  :  when  shock  strength  is  increased,  the  longitudinal  dis¬ 
tance  between  the  shock  wave  and  the  shock  generator  decreases. 
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There  is  therefore  less  "room"  for  downstream  influence. 

Second,  for  a  and  Xp  given,  there  is  no  consistent 
trend  of  the  data  when  changing  unit  Reynolds  numbers.  On  the 
contrary,  for  a  given  unit  Reynolds  number,  there  is  a  consis¬ 
tent  increase  in  upstream  influence  as  Xp  is  increased. 

The  most  useful  observation  of  the  data  being  in  non 
dimensional  form,  a  non  dimensional  plot  of  Lu/6  as  a  function 
of  LSH/6  was  made  (for  a  definition  of  LSH,  see  Fig.  1) 

(Fig.  15).  The  value  of  6  used  in  the  non  d i mens i ona 1 i za ti on 
process  was  the  value  at  the  fin  apex,  as  given  by  the  theory 
of  Chapman  and  Rubesin  (Ref.  24).  This  particular  way  of  presen 
ting  the  data  was  chosen  because  it  proved  to  be  the  suitable 
choice  of  parameters  in  the  turbulent  case  (Ref.  15).  By  dimen¬ 
sional  analysis,  it  can  be  shown  that  with  these  parameters 

chosen,  the  general  formula  for  Lu/6  reads  : 

Lu/<5  =  f  (LSH/6  , Re6  , a)  for  a  fixed  Mach  number. 

In  the  laminar  case  too,  the  data  fit  onto  a  single  band  of 
points.  The  scatter  of  data  is  due  to  uncertainties  in  the 
experimental  values,  which  will  be  discussed  in  chapter  3.4. 
Because  of  this  scatter,  it  is  impossible  to  determine  whether 
there  is  an  influence  of  shock  strength  and  Reynolds  number 
based  on  the  curve  Lu/6  =  f(LSH/fi).  As  an  improvement  of  the 

quality  of  measurements  is  expected  for  the  new  series  of  experi 

ments,  it  is  hoped  that  the  identification  of  the  influence  of 
any  single  parameter  will  be  possible.  A  theoretical  analysis 
similar  to  that  of  Settles  et  al.  (Ref.  26)  will  be  presented 
in  chapter  3.6  in  order  to  predict  the  effect  of  Re6 . 

From  figure  15,  it  can  be  seen  that  nowhere  is  cylin¬ 
drical  symmetry  achieved.  This  does  not  necessarily  mean  that 
the  semi - i nf i ni te  interaction  can  never  reach  cylindrical 
symmetry.  It  just  means  that  in  the  experimental  simulation, 
edge  effects  are  encountered  before  the  hypothetical  appearance 
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of  cylindrical  symmetry.  Yet,  the  quantity  LSH/6  is  already 
large  when  edge  effects  are  susceptible  to  appear  and  it  is 
probable  that  this  type  of  interaction  never  reaches  cylindrical 
symmetry.  Furthermore,  there  is  a  strong  theoretical  argument 
in  this  sense  which  will  be  discussed  in  chapter  3.6. 

Finally,  a  typical  pressure  distribution  is  compared 
to  a  typical  two  dimensional  pressure  distribution  taken  by 
Burgio  (Ref.  21)  at  M  =  2.2  ,  a  =  7.5°,  PQ  =  79  mm  Hg  (which 
under  normal  conditions  of  temperature  gives  Rey  =  1.16  106/m, 

Xp  =  6  cm).  On  figure  16  appears  the  pressure  distribution  by 
Burgio  and  the  pressure  distributions  in  the  three  dimensional 
interaction  for  a  =  8°,  Reu  =  1.31  106/m,  Xp  =  6  m  at  the 
various  lateral  positions.  Major  differences  can  be  observed. 
First,  the  shape  of  the  pressure  distribution  is  radically 
different.  There  are  two  non  exclusive  possible  reasons  for 
this  difference.  The  interaction,  when  viewed  in  a  plane 
perpendicular  to  the  shock  wave  is  that  of  a  boundary  layer 
with  a  "normal"  shock  wave,  not  an  oblique  one.  Also,  the 
introduction  of  three  dimensionality  does  not  simply  reduce 
to  the  existence  of  a  cross  flow  but  modifies  completely  the 
mechanics  of  the  interaction.  Second,  as  lateral  distance  is 
increased,  the  upstream  and  downstream  influences  of  the  inter¬ 
action  become  much  larger  than  that  of  the  two  dimensional 
interaction  yielding  the  same  pressure  ratio.  This  favours  the 
point  of  view  that  the  mechanism  of  the  interaction  is  radically 
different  from  that  of  the  two  dimensional  interaction  :  if  it 
was  similar,  one  would  expect  that,  far  from  the  wedge  apex, 
the  interaction  would  become  quasi  two  dimensional  with  cross 
flow.  As  is  shown  by  a  calculation  using  Lees-Reeves  integral 
method  (Ref.  28)  in  such  a  case,  the  effect  of  cross  flow  is 
small.  Therefore,  the  amount  of  upstream  influence  would  reach 
values  of  the  same  order  as  those  encountered  in  two  dimensional 
interactions  far  from  the  wedge  apex.  This  is  in  contradiction 
with  the  experimental  observation. 
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In  conclusion,  the  major  observations  made  from  the 
pressure  distributions  are  the  following  :  for  all  test  cases, 
there  is  a  considerable  extent  of  the  region  disturbed  by  the 
interaction.  This  extent  is  so  large  that  it  was  not  entirely 
covered  by  the  present  tappings.  It  depends  very  weakly  on 
shock  strength  and  not  measurably  on  Reynolds  number.  The 
pressure  distribution  can  be  divided  into  two  parts  (ahead  and 
downstream  of  the  inviscid  shock  wave)  which  behave  differently. 

3 . 4  Uncertainties  in  the  measurements 

3.4.1  FI  ow^vis^uji^a  t  i  ons^ 

The  flow  visualization  technique  involves  some 
difficulties  due  to  the  nature  of  the  wind  tunnel.  At  Princeton 
(Ref.  15)  a  kerosene  graphite  technique  is  used  with  much  success. 
Unfortunately,  it  was  not  possible  to  use  this  technique  in  the 
S-l  facility  for  the  following  reason.  Since  the  VKI  tunnel 
operates  at  reduced  pressure,  it  needs  a  long  pumping  time 
before  starting  the  flow;  during  this  time  the  kerosene  evapo¬ 
rated.  Therefore,  the  kerosene  had  to  be  replaced  by  a  fluid 
which  would  not  evaporate  during  the  pumping  time.  Oil  was 
selected  because,  as  it  flowed  on  the  model,  it  left  the  graphite 
particles  which  produced  the  visualization.  Nevertheless,  the 
quality  was  much  poorer  than  when  using  the  kerosene  graphite 
technique  because  in  the  latter  case  the  surface  streamlines  are 
much  more  accurately  drawn. 

Furthermore,  due  to  the  inertia  of  the  oil,  it  flowed 
past  the  separation  line  when  injecting  the  model.  At  the 
steady  state,  as  the  test  was  still  being  performed,  the  sepa¬ 
ration  line  could  clearly  be  observed  but  after  shutdown,  there 
was  some  uncertainty  due  to  the  initial  flow  of  the  oil  past 
the  separation  line.  Therefore,  the  method  is  only  qualitative 
and  only  major  trends  can  be  observed. 
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3  .  4  . 2  Surface  p  res^ur  e_measur  ement  s^ 

The  pressure  measurements  have  been  hampered  by 
several  difficulties.  First,  due  to  the  availability  of  the 
data  acquisition  system,  the  number  of  transducers  was  in¬ 
creased  in  order  to  take  the  most  advantage  of  the  micropro¬ 
cessor  abilities.  This  induced  severe  leak  problems  in  the 
tubings,  which  were  not  completely  solved.  Second,  there 
appeared  an  important  zero  shift  of  the  transducers  during 
operation,  so  that  it  was  necessary  to  make  a  zero  correction 
of  the  data.  The  errors  involved  in  the  process  as  well  as 
those  due  to  the  existence  of  small  leaks  are  such  that  the 
pressure  distributions  obtained  at  present  are  only  qualitative. 
They  also  introduce  an  uncertainty  in  the  Mach  number 
determi nati on . 

3 . 5  Laminarity  of  the  interaction 

The  laminarity  of  the  interaction  is  an  important 
issue  because  it  is  essential  that  it  be  ensured  to  allow  com¬ 
parison  with  other  experimental  data  or  with  a  theoretical  pre¬ 
diction.  However,  this  problem  is  much  more  complicated  than 
in  a  two  dimensional  experiment.  In  a  two  dimensional  experi¬ 
ment,  the  occurrence  of  transition  after  the  reattachment 
point  ensured  the  laminar  state  of  the  interaction  because  the 
transition  did  not  exert  any  upstream  influence  is  such  a  case. 
This  could  not  be  checked  experimentally  by  a  direct  method 
but  Ginoux  (Ref.  29)  proposed  an  equivalent  criterion  based 
on  the  variation  of  upstream  influence  with  Reynolds  number. 

For  the  three  dimensional  configuration  under  study,  however, 
transition  can  occur  not  only  longitudinally  but  also  laterally 
in  the  interaction  region.  The  problem  is  then  to  evaluate  the 
extent  of  influence  of  transition  on  the  features  of  the  inter¬ 
action.  This  problem  was  addressed  by  Korkegi  (Ref.  28)  in  the 
case  of  a  supersonic  corner.  It  was  noticed  that  a  characteris¬ 
tic  bend  (Fig.  7)  occurs  in  the  separation  line  at  the  location 
where  transition  takes  place.  This  was  never  encountered  in 


17 


the  present  test  series,  not  even  for  a  =  10  or  153.  As  a 
matter  of  fact,  the  experimental  conditions  were  chosen  such 
that  an  equivalent  two  dimensional  interaction  would  be  in  a 
fully  laminar  state.  That  observation  thus  strongly  supports 
the  laminar  character  of  the  tests  performed,  since  there  is 
no  apparent  effect  on  the  surface  flow  visualization.  In 
addition,  the  very  large  extent  of  separated  flow  supports  the 
laminar  nature  of  the  test  and  the  weak  upstream  influence  of 
transition.  Nevertheless,  as  long  as  there  is  no  theoretical 
prediction  of  the  upstream  extent  of  influence  of  transition, 
it  is  difficult  to  assess  exactly  which  part  of  the  flow  can 
actually  be  considered  of  the  laminar  type. 

3 . 6  Discussion  of  the  results 

In  this  section,  the  question  of  the  effect  of 
Reynolds  number  and  the  question  of  cylindrical  symmetry  will 
be  discussed  as  announced  in  chapter  3.3.  In  addition,  a  ten¬ 
tative  explanation  of  the  differences  between  two  dimensional 
and  three  dimensional  interactions  will  be  proposed. 

The  Reynolds  number  effect  on  three  dimensional  shock/ 
turbulent  boundary  layer  interactions  was  studied  in  two  papers 
by  Settles  et  al .  (Ref.  16)  and  Dolling  (Ref.  9),  using  dimen¬ 
sional  analysis.  Starting  from  the  two  dimensional  formula  for 
upstream  influence  as  a  function  of  various  parameters 

0.23a  -1/3 

Lu/S  =  0 . 9e  Res  for  M  =  3  (1) 

they  looked  for  a  three  dimensional  generalization.  The 
general  formula  for  upstream  influence  in  the  skewed  shock 
wave  boundary  layer  interaction  being 

Lu/«  =  f  (LSH/6  ,  Re6  ,a,M)  (2) 

they  assumed  that  the  dependence  on  Reynolds  number  in  the 
three  dimensional  case  would  exhibit  the  same  form  as  in  the  two 
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dimensional  case,  namely  they  restricted  the  most  general 
expression  to  a  single  one 


Lu  , 

<5 


a 

Reo 


LSH 
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b  1 
/  Re 5  ,  a  ,M  j 
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When  compared  to  the  experimental  data  this  expression  was  able 
to  correlate  the  experimental  data  and  a  and  b  had  the  common 
value  of  1/3  which  was  already  used  for  two  dimensional  scaling 
Moreover,  it  was  found  that  the  upstream  influence  was  not  a 
function  of  the  two  independent  variables  M  and  a  but  only  of 
a  combination  of  the  two,  namely  the  normal  Mach  number 
Mn  =  M s i n 3  (for  M  and  n  fixed,  3  is  known).  If  an  analogous 
development  is  performed  in  the  laminar  case,  the  two  dimen¬ 
sional  relation  being  : 


Lu 


3/2 

SRe<5 


f(M,a) 


(4) 


the  three  dimensional  relation  will  read  : 


6  R  e  * 

This  relation  remains  to  be  checked  experimentally. 

The  cylindrical  symmetry  of  the  interaction  is  also 
a  problem  which  has  been  posed  a  long  time  ago  by  McCabe 
(Ref.  1)  in  the  turbulent  case.  Examining  results  in  a  similar 
geometry,  Nebbeling  (Ref.  30)  concluded  that  the  interaction 
could  be  divided  in  several  regions  (Fig.  18).  From  that 
figure,  it  can  be  stated  that,  even  if  the  separation  line 
runs  parallel  to  the  inviscid  shock  wave,  the  flow  may  never  be 
purely  of  the  cy 1 i ndr i ca 1 1 y  symmetric  type  since  the  distance 
between  the  shock  wave  and  the  wedge  keeps  increasing  as  one 
progresses  laterally  (see  "conical  region"  behind  "two 
dimensional  region"  in  Fig.  18). 
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As  discussed  earlier,  important  differences  appear 
between  the  behaviour  of  two  dimensional  and  three  dimensional 
interactions.  In  particular,  the  appearance  of  separation 
and  the  amount  of  upstream  influence  are  two  important  features 
different  in  both  types  of  interactions.  This  is  due  to  the 
fundamental  difference  in  nature  between  two  dimensional  and 
three  dimensional  separation.  Recently,  the  issue  of  three 
dimensional  separation,  which  has  been  controversial  since  the 
pioneering  papers  of  Maskell  (Ref.  31)  and  Ei chel brenner  et  al. 
(Refs.  32,33),  has  been  clarified  by  a  significant  contribution 
from  Peake  and  Tobak  (Refs.  17,22,34).  From  continuity  argu¬ 
ments  and  topological  considerations,  the  nature  of  separation 
in  flows  over  simple  configurations  was  established  and  flow 
models  presented.  A  major  difference  which  exists  between  two 
dimensional  and  three  dimensional  separated  flows  clearly 
appears  on  figure  19  where  are  shown  possible  flow  patterns  in 
the  plane  of  symmetry  ahead  of  a  three  dimensional  obstacle  and 
two  dimensional  flow  pattern  ahead  of  a  step.  The  major  dif¬ 
ference  can  be  stated  as  follows  :  in  two  dimensional  flow,  it 
is  impossible  to  draw  a  path  from  the  separated  region  to  the 
upstream  infinity  without  crossing  the  dividing  line  whereas 
in  three  dimensional  flow,  the  "separated  region"  is  continuously 
fed  by  fluid  coming  from  upstream.  A  most  spectacular  example 
of  this  is  the  flow  over  a  delta  wing.  To  state  the  same  thing 
in  another  way  :  in  two  dimensional  flow,  a  separated  region 
is  really  separated  from  upstream  infinity  by  the  dividing  line, 
whereas  in  three  dimensional  flow,  this  is  untrue  :  as  pointed 
out  by  Korkegi  (Ref.  4),  there  is  a  certain  amount  of  mass 
transfer  toward  the  separated  region.  Tne  importance  of  this 
mass  transfer  is  of  greatest  importance  in  the  nose  region  of  a 
blunt  fin  induced  shock  wave  boundary  layer  interaction  as 
mentioned  in  reference  35.  Therefore,  the  word  separation  does 
not  carry  the  same  sense  in  two  dimensional  and  three  dimen¬ 
sional  cases.  A  consequence  of  the  mass  transfer  is  that,  if 
the  reattaching  streamline  comes  from  outside  the  incoming 
boundary  or  from  high  energy  layer  in  the  boundary  layer,  the 
boundary  layer  which  develops  behind  the  separation  line  is  a 
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"new"  boundary  layer,  independent  of  the  incoming  boundary 
layer  character i sties .  This  is  the  reason  for  the  very  high 
heat  transfer  rates  encountered  in  three  dimensional  hypersonic 
separated  flows. 

Yet,  there  are  examples  of  three  dimensional  flows  in 
which  separation  exhibits  features  very  close  to  that  of  two 
dimensional  separated  flows,  an  example  of  which  is  shown  in 
figure  20,  together  with  the  equivalent  two  dimensional  configu¬ 
ration.  In  this  case,  the  mass  transfer  towards  the  separated 
region  is  so  small  that  the  flow  may  be  considered  of  the  two 
dimensional  type.  This  is  of  course  the  basis  for  carrying  out 
two  dimensional  experiments. 

The  previous  discussion  paves  the  way  for  establishing 
a  criterion  for  the  classification  of  the  various  three  dimen¬ 
sional  s hoc k/bou ndary  layer  interactions.  It  is  obvious  that 
the  mechanism  of  three  dimensional  separated  flows  is  radically 
different  in  the  cases  of  large  or  very  small  mass  adjunction 
in  the  separated  region.  Therefore,  a  qualitative  estimation 
of  the  mass  adjunction  would  provide  a  suitable  criterion  for 
distinguishing  the  various  three  dimensional  separated  flows. 
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4 •  theoretical  investigations 

As  mentioned  in  the  introduction,  a  theoretical 
investigation  was  performed.  Two  approaches  were  studied  in 
more  detail  :  the  newly  proposed  corrected  viscosity  method 
for  solving  the  full  Navier  Stokes  equations  and  a  possible 
extension  of  the  Lees-Reeves  integral  method. 

For  solving  these  kinds  of  flow  problems  where  there 
is  a  strong  coupling  between  viscous  and  inviscid  flow  there 
are  essentially  two  approaches.  The  first  one  consists  in 
solving  the  full  Nav i er-Stokes  equations  while  the  second 
consists  of  performing  a  coupled  computation  of  an  inviscid 
external  flow  field  with  a  boundary  layer. 

The  first  type  of  method  uses  finite  difference  or 
finite  element  techniques  for  discretizing  the  Navi er-Stokes 
equations.  This  renders  the  methods  very  general,  so  that 
after  a  major  initial  effort,  many  flow  configurations  could  be 
studied.  To  this  category  belong  the  methods  of  McCormack 
(Ref.  12),  Beam  and  Warming  (Ref.  13),  and  the  implicit  cor¬ 
rected  viscosity  method  recently  proposed  by  a  VKI  faculty 
member.  Prof.  Essers  (Ref.  14)  (the  latter  method  is  studied 
in  further  detail  in  chapter  4.1J.  Their  main  disadvantage 
is  that  they  require  considerable  storage  space  and  long 
computing  times. 

Among  methods  using  a  coupling  approach  (which  are 
discussed  in  detail  by  Le  Balleur,  Ref.  36),  the  integral 
method  in  general  and  the  integral  method  of  Lees-Reeves 
(Ref.  11)  in  particular  offer  the  advantage  of  a  reduced  storage 
requirement  or  computing  time,  with  the  disadvantage  of  a  small 
flexibility.  Since  the  VKI  has  considerable  experience  with 
these  methods,  the  possibility  of  adapting  the  work  of  Leblanc 
(Ref.  28)  to  the  present  flow  configuration  was  studied 
(Ch  a  p t .  4.2). 
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4 . 1  Implicit  corrected  viscosity  method  for  solving 
the  compressible  Na v i er-S tokes  equations 

Presently,  the  two  most  effective  methods  for  solving 
the  full  compressible  Nav i er- Stokes  equations  are  McCormack's 
(Ref.  12)  and  Beam  and  Warming's  (Ref.  13).  The  main  advantage 
of  McCormack's  method  is  that  it  needs  the  solution  of  block 
bidiagonal  systems  whereas  block  tridiagonal  systems  should  be 
solved  in  Beam  and  Warming's  method.  On  the  other  hand.  Beam 
and  Warming's  method  has  the  big  advantage  that  the  precision 
of  any  steady  state  solution  is  independent  of  the  time  step 
used  in  the  marching  process. 

Recently,  new  methods  based  on  the  principle  of 
implicit  corrected  viscosity  were  proposed  (Ref.  14).  Their 
principle  is  as  follows  -  suppose  that  the  steady  state  equa¬ 
tions  to  be  solved  are  the  following  : 

D.  (S)  =  0  i  =  1,  ...  M  (1/ 

D.  being  space  differential  equations  applied  to  the  unknown 
vector  S.  They  correspond  to  the  steady  part  of  compressible 
Navier-Stokes  equations  or  a  combination  thereof.  The  unsteady 
equations  solved  in  all  classical  time  dependent  techniques  can 
be  written  as  : 


i  ■*  ■* 

—  =  MS)  +  xv. (S) 
at 


(2) 


where  X  is  a  small  positive  constant  and  the  V ^ ' s  are  space 
elliptic  aritificial  viscosity  operators  possibly  used  to  avoid 
non  linear  instabilities,  when  hyperbolic  equations  are  to  be 
solved.  As  the  artificial  viscosity  introduces  errors  to  the 
steady  solution,  X  should  obviously  be  kept  small.  This  limits 
the  speed  of  convergence  of  the  method.  To  significantly 
improve  the  rate  of  convergence,  Couston,  McDonald  and 
Smolderen  (Ref.  37)  proposed  a  so-called  corrected  viscosity 


approach.  In  that  technique,  the  following  equations  are 
solved  : 

3  w  - 

—  =  Di(s)  +  x|v.(s)-v.(s*) ;  (3) 

Yt  c-  - 

-►  ^  -> 

S  being  the  vector  S  computed  at  some  previous  iteration. 

At  steady  state,  S  =  S*  and  therefore,  the  accuracy  of  the 
solution  does  not  depend  on  the  quantity  \ .  On  the  other  hand, 
the  introduction  of  the  correction  term  does  not  alter  the  dam¬ 
ping  qualities  of  equation  (2).  The  method  is  therefore  expec¬ 
ted  to  yield  faster  convergence  speeds  than  existing  methods. 
This  has  proved  to  be  the  case  for  explicit  corrected  viscosity 
techniques  which  have  been  successfully  used  for  various  tran¬ 
sonic  problems  including  aerofoil  and  cascade  flows. 

What  is  proposed  here  is  to  use  implicit  corrected 
viscosity  schemes.  Essers  has  proved  (Ref.  14)  that  the  use  of 
a  one  step  fully  implicit  scheme  is  not  interesting.  To  over¬ 
come  this  difficulty,  it  is  suggested  to  perform  a  splitting  of 
equations  (3).  In  that  technique,  a  series  of  similar  cycles 
are  computed  successively  at  various  time  levels.  Each  cycle 
consists  of  an  explicit  physical  part  followed  by  an  implicit 
corrected  viscosity  step.  This  has  the  further  advantege  that, 
if  complicated  turbulence  models  have  to  be  introduced  in  the 
physical  equations,  they  would  be  introduced  explicitly  and 
would  not  require  the  calculation  of  complicated  jacobians  as 
required  by  the  methods  of  Beam  and  Warming  or  McCormack.  The 
explicit  schemes  could  be  composed  of  several  steps  but  should 
be  dissipative  for  all  wavelengths.  The  corrected  viscosity 

a  £  +  1 

steps  consist  of  correcting  the  explicitly  predicted  value  S 
by  solving  the  following  equation  : 


This  could  be  done  by  overrelaxation  methods  or  more  efficiently 
by  multigrid  techniques. 
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This  method  was  studied  and  tested  over  a  simple  one 
dimensional  test  case  (Ref.  38).  It  proved  to  be  efficient. 

The  extension  to  two  dimensional  cases  was  studied  theoretically 
using  a  Von  Neumann  analysis.  The  results  are  presented  in 
detail  in  reference  34.  Let  us  here  briefly  state  the  major 
findings.  Suppose  that  the  equation  to  be  solved  reads  : 


3JJ  +a  a_u  a  2u 

at  ax  ax 
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it  is  a  one  dimensional  equation  involving  both 


conduction 
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3  X 


and  diffusion 
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It  thus  figures  the 


the  features  of  fluid  flow. 


Among  explicit  schemes  which  can  be  used  to  compute 
the  physical  step,  the  Lax  Wendrof f  2-step  scheme  exhibits  a 
major  advantage.  If  it  is  used  to  solve  equation  (5)  the 
scheme  wi 1 1  read  : 
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Corrected  viscosity  : 
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It  has  been  proved  (Ref.  33)  that  the  faster  convergence  speed 
is  obtained  for  uQ  =  u\  which  means  that  the  viscous  term 
should  be  corrected  by  the  previous  time  step  value. 

In  such  a  case,  the  latter  two  steps  can  be  combined 
and  written  in  A  form. 

1  +  2  i 
tsu  l  =  u  -  u 


z+1 


z  +  1 


All 


2  a  A  t 


3  u 


3xm 


+  2  e  At 


3^u 


3x  2m 


+  2y  At 


.u 


3x~m 


This  proves  one  advantage  of  the  scheme  (shared  by  that  of 
Beam  &  Warming)  :  the  accuracy  of  the  converged  solution  (time 
level  i+1)  does  not  depend  on  the  time  step  used. 

When  a  stability  analysis  of  the  scheme  is  performed 
(Ref.  39)  it  is  shown  that  with  suitable  discretization  of  the 
spatial  derivatives,  the  stability  and  convergence  speed  of 
the  scheme  depends  on  two  parameters 

0  =  X  =  — 

u  a 

where  j  is  the  pulsation  of  the  mode  under  study.  For  each  \ 
there  is  an  optimum  value  of  e  (-1)  which  maximizes  the  conver¬ 
gence  speed.  On  the  contrary,  it  appears  that  there  is  no  limit 
on  a.  For  increasing  \,  one  should  get  increasing  convergence 
speed.  This  will  probably  be  limited  by  boundary  conditions 
or  non  linear  instabilities. 


Yet,  a  convergence  speed  of  order  unity  may  be  obtained 
using  such  a  method.  The  extension  to  two  dimensional  problems 
can  easily  be  done.  The  qualities  of  the  method  are  conserved. 
This  very  encouraging  result  ranks  the  method  among  the  most 
promising  ones. 
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4 . 2  Extension  of  Lees-Reeves  integral  method 

Following  the  works  of  Leblanc  (Ref.  23)  who  introduced 
the  effect  of  cross  flow  in  an  axisymmetric  interaction  between 
a  shock  wave  and  a  laminar  boundary  layer,  it  was  envisaged  to 
extend  this  method  to  compute  a  cy 1 i ndr i ca 1 1 y  symmetric  normal 
shock/laminar  boundary  layer  interaction  with  cross  flow  which 
would  represent  the  cyl i ndrical ly  symmetric  part  of  a  skewed 
s hoc k/ 1  am i na r  boundary  layer  interaction. 

However,  there  is  a  major  problem  because,  the  inviscid 
flow  behind  the  shock  being  subsonic  (in  the  normal  direction), 
the  coupling  used  by  Lees  and  Reeves  (Prandtl  Meyer  relation) 
could  not  be  used.  A  subsonic  coupling  would  by  no  means  be 
impossible  but  would  be  rather  involved.  Furthermore,  the 
experimental  results  show  no  region  of  cylindrical  symmetry  for 
thi s  i nteraction . 

On  the  other  hand,  the  method  appears  well  suited 
for  the  computation  of  interactions  at  swept  compression  corners 
which  is  known  to  exhibit  regions  of  cylindrical  symmetry  in 
the  turbulence  case  (Ref.  26).  This  configuration  seems 
analogous  to  that  studied  by  Leblanc  (Ref.  28)  but  important 
differences  exist: in  this  case,  the  total  enthalpy  is  constant 
throughtout  the  boundary  layer  (for  Pr  =  1).  Also,  if  sweep 
is  to  be  introduced  for  a  fixed  Mach  number,  the  normal  Mach 
number  will  vary,  which  is  not  the  case  in  the  work  of  Leblanc. 
This  is  the  reason  why  the  analysis  of  the  method  was  persued. 

A  preliminary  study  was  carried  out  whose  results 
are  presented  in  reference  40.  A  summary  of  the  analysis  is 
given  in  Appendix  I.  This  shows  that  the  effect  of  sweep  is 
anal ogous  to  the  effect  of  heating.  It  is  therefore  expected 
that  it  will  have  a  similar  effect  as,  for  instance,  increasing 
the  upstream  influence. 
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A  computation  of  the  weak  interaction  region  which 
develops  at  the  leading  edge  of  a  plate  in  supersonic  flow  was 
carried  out  in  order  to  verify  this  prediction.  The  results 
(Fig.  21)  show  that  it  is  verified.  This  partly  explains  the 
unfavourable  effect  of  sweep  on  the  extent  of  disturbed 
regions  in  shock/boundary  layer  interaction  problems. 
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5 .  POSSIBILITY  OF  PERFORM  I N G  LDV  MEASURE M ENTS 
IN  THE  SKEWED  SHOCK  WAVE/ LAMINAR  BOUNDARY 
LAYER  INTERACTION 


As  announced  in  the  original  proposal,  the  possibility 
of  performing  LDV  measurements  in  this  interaction  was  evaluated. 
It  appeared  that  the  measurements  would  present  considerable  dif¬ 
ficulty  which  would  prevent  obtaining  results  in  the  near  future. 
Yet,  a  progressive  program  aiming  at  performing  velocity  measure¬ 
ments  in  this  interaction  can  be  constructed.  This  is  the  object 
of  a  new  research  proposal  to  AFOSR . 
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6 .  CONCLUSIONS 

An  experimental  investigation  of  a  three  dimensional 
skewed  shock  wave/1  ami nar  boundary  layer  interaction  has  been 
performed.  Surface  flow  visualizations  and  surface  pressure 
measurements  have  been  carried  out  over  a  variety  of  test  condi¬ 
tions.  This  provided  a  parametric  study  of  the  interaction. 
Simultaneously,  two  theoretical  approaches  were  evaluated  :  an 
extension  of  Lees  -  Reeves - K1 i neb  erg  integral  method  and  new 
implicit  corrected  viscosity  methods.  The  conclusions  are  : 

The  upstream  influence  is  much  larger  in  this  three  dimen¬ 
sional  interaction  than  in  a  two  dimensional  interaction 
having  the  same  overall  pressure  rise.  So  is  the  downstream 
influence.  This  had  the  consequence  that  the  pressure  tappings 
did  not  cover  the  entire  interaction  region. 

Extended  separation  occurs  even  for  the  smallest  wedge 
incidence,  it  is  for  a  pressure  ratio  of  1.27. 

A  criterion  is  proposed  to  characterize  three  dimensional 
interactions  and  to  distinguish  between  quasi  two  dimensional 
interactions  and  fully  three  dimensional  interactions,  the 
mechanisms  of  each  type  being  different. 

The  integral  method  of  Lees-Reeves-Kl i neberg  is  not  suited 
for  the  analysis  of  the  skewed  shock  wave/laminar  boundary  layer 
interaction.  On  the  contrary,  it  is  well  suited  for  the 
analysis  of  the  flow  at  swept  compression  corners.  It  indicates 
that  sweep  tends  to  increase  upstream  influence. 

The  newly  proposed  implicit  corrected  viscosity  method  was 
studied.  The  analysis  shows  that  by  suitably  choosing  its 
parameters,  convergence  speeds  or  order  unity  can  be  achieved. 
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APPENDIX  I  -  LEES  REEVES  INTEGRAL  METHOD  FOR 
SHOCK/BOUNDARY  LAYER  INTE RACE  I  ON 
WITH  CROSS  FLOW 

Under  the  assumption  of  cylindrical  symmetry  and 
adiabatic  conditions,  the  boundary  layer  equations  read  : 

Continuity  : 

—  (PU)  +  ---  (;.V)  =  0 
3  x  3y 

x  - momentum  : 


3  u  3  u 

p  u  —  +r  v  — 


3x  3y  3x 

z-momentum  : 


.  i£  +  ^  [, 

3  x  3  y  ^  3  y  j 


3w  aw  a  f  awl 

PU  -  +pv  -  =  -  p  - I 


ax  ay  3y  {  3 y] 

One  should  notice  the  analogy  between  this  latter  equation  and 
the  energy  equation  (Pr  =  1). 


3  H 

3  H 

3 

3  H  1 

pa  —  + 

pv  — 

“  —  I 

3x 

ay 

ay 

ay  j 

This,  under  adiabatic  conditions  has  the  trivial  solution 
H  =  He  =  H  w 


Introducting  the  following  notation 
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JW  = 


=  constant 


V  2  H  e 


•  w  - 
1-  -w: 


K  =  ^ 
we 


and  integrating  the  equations  across  the  boundary  layer,  with 
the  help  of  the  II 1 i ngsworth-Stewartson  tr a nsf orma t i on ,  one 
gets  the  following  integral  equations. 


F  -Ji  +5*  (if  ii  +  if  is;  +f o*  i_ 

dx  1  [sa  dx  3 C  dxj  1  dx 


— 1  +£*  if  ii  +£*  (1+2H+UEPS)  -d—  r<nMe  j  =  6 
dx  ds  dx  1  rtu  t  ai 


d 

dy 


mMe  ]  =  SC 

Uj  h 


=  6C. 


d6  .  ,  ,  , 

j  — 1  +5*  if  ii  +5*  (3J-2UTAU|  —  ftnMe  ) 
1  J-  J'  1  dx  ^  “J 


dx 


da  dx 


'■  d  f  tnMe  1  =  3CU  ^ 


M 

h 

Me 

a 

Re  &*. 
a  1 

M 

P 

M 

cr>  'x 

Re  f* 

a  1 

M 

ex-  a 

R 

di, 

Tw  — 1  +<5* 
dx  1 


3  Tw  da  a  Tw  dc 
3a  ax  3c  dx 


h 

d  c .  ^  1  '®cx  Qw 


Me  Re  6T 

a  a  i 

M 


+  5.Tw  —  (tnMe)  =  gC  - - 

dy  "  Me  Re  &*. 

a  a  1 


Me^  and  Re^  denote  Mach  and  Reynolds  number  based  on  the  exter¬ 
nal  longitudinal  velocity  ue  and  u^.  The  quantities  appearing 
in  the  equations  are  presented  in  Table  2.  The  integral  quan¬ 
tities  can  be  expressed  as  integrals  of  non  dimensional 
longitudinal  and  transverse  velocity  profiles,  for  example 


EPS  -  a  (1-K2 )  dy 


n  -|*i 

|  (l-f')dyj 

i  o  .! 


where 
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The  identification  parameters  are  : 


a  =  rSf ' ' (0) 


for  attached  flow  C  =  K'(0) 


nf  '  =0 

n  5 


for  separated  flow 


The  velocity  profiles  (f'  and  K)  are  chosen  to  be  the  similar 
solutions  of  the  boundary  layer  equations.  They  are  thus 
the  solutions  of  the  Falkner  Skan  equations  : 


\f ' '  '  +  f f ' ’ +  p  [  1  - f ' 2  +  u(l-K2)  l  =  0 

< 

K' ' +f K '  =  0 

with  the  boundary  conditions 


n  =  0 


f  =  f  ’  =  0 


K  =  0 


n 


no 


i  im 

f'  =  1 

r|-»a> 


n  m 

K  =  1 

n-*-“ 


It  can  be  noticed  that  the  integral  quantities  do  not  depend 
explicitly  on  the  parameter  oj  but  only  through  the  dependence 
on  u  of  the  functions  f  and  K.  This  dependence  can  be  expected 
to  be  small  and  a  unique  representation  for  the  integral  quan¬ 
tities  as  a  function  of  the  identification  parameters  will  be 
posible.  This  has  been  verified. 


The  coupling  is  given  by  a  Prandtl  Meyer  relation  in 
which  all  quantities  are  considered  in  the  longitudinal  plane. 


o 

e 


- 1 


tg 


ve 


u 


e 


o0  + 


N 


rj£ r 

Y  -  1 


(M^-l)-tg'1 


a  :  P  =  11000  Pa 
b  :  P  =  22000  Pa 

All  distances  in  millimeters 


TABLE  1  -  UPSTREAM  INFLUENCE 
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TABLE  2 


FIG  3 -MODEL  INJECTED  IN  THE  TUNNEL  FOR  SURFACE 
FLOW  VISUALIZATION  PHOTOGRAPHED  AFTER 


LOW  VISUALIZ 


MACH  CONE  FROM  POINT  A 


NOTATIONS  AND  RELATIONS 


a  :  WEDGE  INCIDENCE 
(3  :  SHOCK  WAVE  ANGLE. 

MACH  NUMBER  IN  THE  FREE 


Ml  : 


STREAM 


MACH  NUMBER  BEHIND 
OBLIQUE  SHOCK 


AO  =  r 

OA'  r  rtg  (  P  -O  ) 

OB  =  rtg  M- 2  WITH  M-2  =  Sin  _1(  —  ) 

AA’=  — _ 

cos ( p  -  oc) 

—  coS4,oa;=  a<Mi> 

OB  tg  |12 

A'B  =  r  tg  M-2  Sin  © 

tg V  =  =  tg  M-2  Sin9cos(P- 

A  A' 


FIG.  B- EXPERIMENTAL  CONFIGURATION  WITH  SHOCKWAVE 


FIG.  12  -  STATIC  PRESSURE  DISTRIBUTIONS, Y  =  2.25  INCH  ,  MODEL  2  {from  ref.  6) 


from  ref.  6  ) 


AM  INFLUENCE  AS  A  FUNCTION  OF  DISTANCE  FROM  THE  PIN 


IONS  IN  A  3  D  SKEWED  SHOCK  WAVE  TURBULENT  BOUNDARY 


S'  :  SEMI -SADDLE  POINT,  N  .FOCUS  OR  SPIRAL  NODE 
(NOTATIONS  TAKEN  FROM  Ref.  17  I 

POSSIBLE  FLOW  PATTERN  IN  PLANE  OF  SYMMETRY 
AHEAD  OF  A  CYLINDRICAL  PROTUBERANCE 
(AFTER  FLOW  VISUALIZATIONS  IN  Ref.  35  ) 


S' 


C  .CENTRE 

FLOW  PATTERN  AHEAD  OF  A  2  D  STEP 


FIG.  19-3D  &  2D  SEPARATED  PATTERNS  . 


b)  FLOW  PATTERN  IN  THE  CUT  IN  THE  PLANE  OF  SYMMETRY 


FIG.  20  - 


(QUASI  2  D)  SEPARATED  FLOW  PATTERN 


REX 

FIG.  21  -  EFFECT  OF  SWEEP  ON  THE  WEAK  INTERACTION 
REGION 
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